Rogue waves in the Davey- Stewart son equation 
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General rogue waves in the Davey-Stewartson-I equation are derived by the bilinear method. It 
is shown that the simplest (fundamental) rogue waves are line rogue waves which arise from the 
constant background with a line profile and then disappear into the constant background again. It 
is also shown that multi-rogue waves describe the interaction of several fundamental rogue waves. 
These multi-rogue waves also arise from the constant background and then decay back to it, but in 
the intermediate times, interesting curvy wave patterns appear. However, higher-order rogue waves 
are found to show more interesting features. Specifically, only part of the wave structure in the 
higher-order rogue waves rises from the constant background and then retreats back to it, and this 
transient wave exhibits novel patterns such as parabolas. But the other part of the wave structure 
comes from the far distance as a localized lump, which decelerates to the near field and interacts 
with the transient rogue wave, and is then reflected back and accelerates to the large distance again. 
These rogue-wave solutions have interesting implications for two-dimensional surface water waves 
in the ocean. 



I. INTRODUCTION 

Rogue waves are large and spontaneous ocean surface 
waves that occur in the sea and are a threat even to 
large ships and ocean liners pQ. Recently, an optical 
analogue of rogue waves — optical rogue waves, was 
also observed in optical fibres [2j [3]. A growing con- 
sensus is that both oceanic and optical rogue waves ap- 
pear as a result of modulation instability of monochro- 
matic nonlinear waves. Mathematically, the first and 
simplest rogue-wave solution was reported in the non- 
linear Schrodinger (NLS) equation by Peregrine [4 . This 
solution approaches a non-zero constant background as 
time goes to ±oo but develops a localized hump with peak 
amplitude three times the constant background in the 
intermediate times. Recently, higher-order rogue waves 
in the NLS equation were reported in many articles [5]- 
[T2] . It was shown that these higher-order waves could 
reach higher peak amplitudes or exhibit multiple inten- 
sity peaks at different spatial locations and times. In 
addition to the NLS equation, rogue waves in some other 
wave equations (such as the Hirota equation) have also 
been explored [13]. Rogue waves are intimately related 
to homoclinic solutions which approach a constant back- 
ground as time goes to ±oo but develop spatially-periodic 
wave patterns in the intermediate times [T4UT7 j . Indeed, 
rogue waves can be obtained from homoclinic solutions 
when the spatial period of homoclinic solutions goes to 
infinity M CD3 EE] • 

Rogue waves which have been studied so far are mostly 
one-dimensional. But ocean surface waves are always 
two-dimensional. Thus a natural question is to inves- 
tigate rogue waves in two-dimensional model equations. 
It is well known that the evolution of a two-dimensional 
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wave-packet on water of finite depth is governed by 
the Benney-Roskes-Davey-Stewartson equation [T9H2T] . 
This equation is often just called the Davey- Stewart son 
(DS) equation in the literature. The DS equation is di- 
vided into two types, DSI and DSII equations, depending 
on the signs of its coefficients. The simplest (one-mode) 
homoclinic solution to the DS equation was derived in 
[22] . Taking the spatial period of this homoclinic solution 
to go to infinity, the simplest (fundamental) rogue-wave 
solution was also obtained there. But more general rogue 
waves in the DS equation are still unknown. 

In this paper, general rogue waves in the DS-I equa- 
tion are derived. These solutions are obtained by the 
bilinear method and expressed in terms of determinants. 
It is shown that the simplest (fundamental) rogue waves 
are line rogue waves which arise from the constant back- 
ground with a line profile and then disappear into the 
constant background again (this simplest rogue wave 
agrees with that reported in [22]). It is also shown 
that non-fundamental rogue waves contain different types 
such as the multi-rogue waves and higher-order rogue 
waves. The multi-rogue waves describe the interaction 
of several fundamental rogue waves. These multi-rogue 
waves arise from the constant background and then de- 
cay back to it, but in the intermediate times, interesting 
curvy wave patterns appear. Higher-order rogue waves, 
on the other hand, exhibit certain features which are very 
novel. Specifically, only parts of the wave structures in 
these higher-order rogue waves rise from the constant 
background and then retreat back to it, exhibiting un- 
usual transient wave patterns (such as parabola shapes) 
in the intermediate times. But the other parts of the 
waves come from the far distance as localized lumps, 
which interact with the transient rogue waves in the near 
field and then are reflected back to the large distance 
again. Since the DS equation is a well known mathe- 
matical model for two-dimensional surface water waves, 
these general rogue-wave solutions can have interesting 
implications for the study of rogue waves in the ocean. 
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II. RATIONAL SOLUTIONS IN THE 
DAVEY- STEWARTS ON- 1 EQUATION 

The Davey-Stewartson-I (DSI) equation is given by 



iA t = A xx + A yy + (e\A\ 2 - 2Q)A, 
Qxx Qyy ^(l^-l )xxi 



(2.1) 



where e = 1 or — 1. It is noted that under the variable 
transformation Q — >• Q + e|A| 2 , x ^ y and e — » — e, 
this equation is invariant, thus we can fix the sign of e 
without loss of generality. However the transformation 
Q —> Q + e\A\ 2 changes the boundary condition of Q in 
general, thus we keep e in our analysis. Equation (2.1) is 
transformed into the bilinear form, 



(Dl + Dl-iD t )g.f = 0, 
(Dl-Dl)f.f = 2e(f*-\g\*), 

through the variable transformation, 

A = y/2^, Q = e-(2Iog/) xx , 



(2.2) 



(2.3) 



where / is a real variable and g is a complex one. 

Rogue waves are rational solutions (under certain pa- 
rameter restrictions). Thus we first present the general 
rational solutions to the DSI equation in the following 
theorem. The proof of this theorem is given in the ap- 
pendix. 



Theorem 1 The DSI equation (2.1) has rational solu- 



tions (2.3) with / and g given by N x N determinants 
f = r , g = n, (2.4) 

where r n = deti<^j<iv ( m l^)' anc ^ ^ e matrix elements 

are given by either 

(a) 



i n) = E 

k=0 

rij 

E 

1=0 



c ik (Pid Pt + £ + n) n '- 



Tin— I 



1 



Pi + Pj 



(2.5) 



£ = Vi ~f^ + Pl+ : P7 \ + P ' + r^ 2 t, (2.6) 



or 
(b) 



(n) 



ni+rij 



(2.8) 



where S n (x) is the elementary Schur polynomial defined 
via the generating function 



J2S n (x)\ n = exp [J2x k \ k 



n=0 



\k=l 



for x = (xi,x 2 , • • • ), 

£ (n) (p) = [£i(p)+™, 6(p), 



ffc(p) + fain, 



Sij is the Kronecker delta notation (which is equal to 1 
when i = j and zero otherwise), and 



&(p) = ^ 



p + e(-l) fc /p p-e(-l) fc /p 



2V-(-2) fc /p 2 
2a/^T 



In (a) and (b), p^, and are arbitrary complex con- 
stants, and rii is an arbitrary positive integer. These two 
expressions in (a) and (b) would yield identical solutions 
if the constants and in them are related by 

Cik = Cikfai -k)\, dji = dji(rtj -1)1. 

Thus in the later text we use the expression in (a). By 
a scaling of / and g, we can normalize Qo = 1 without 
loss of generality, thus hereafter we set Qo = 1. We will 
also call the above solution as the TV-rational solution of 
order (ni, 77,2, • • • , tin). We comment that a more explicit 
expression for mf^ similar to (2.6) in [12] can also be 
obtained, but since that expression is a bit complicated, 
we omit it in this paper. 

The simplest rational solution, namely 1-rational solu- 
tion of 1st order, is given by taking N = 1 and n\ = 1, 
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Pi +Pl 
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Pi 



Pi +Pl 
_ Pi 
Pi + P1 
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(Pl+Pl) 2 
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where 

, _ pi - epi 1 pi + cpr 1 .. , Pj+Pl 2 + 

Si- 2 X+ 2 y+ i ' 

and pi, en are arbitrary complex constants. This solu- 
tion can be rewritten as 

1 -(ee+A), ^ = ^^[(e+i)(e-i) + A], 



pi + pi 

where 



Pi + Pi 



and 



£ = ax + by + ujt + 0, A = PiPi/(pi + Pi) 2 , 



a = (pi - e^ 1 )/2, 6 = (pi + e^ 1 )/2, 



cj = (Pi +Pi 2 )/i, = en -pi/(pi +Pi). 
If we separate the real and imaginary parts of a, 6, uo and 



From the above analysis, we see that rational solutions 



(2.4) to the DSI equation become rogue waves when the 
parameters Pi are real- valued (this fact holds for N = 1 
as well as for higher N integers). In the next section, we 
will examine these rogue waves in more detail. 

It is noted from the above explicit solution formulae 
that the parameter en causes a shift of the origin in the 
(x,y,t) space. Thus en can be set to zero by a shift of 
the (x,y,t) axes. This fact also holds fo r TV = 1 as well 
as for higher N integers in the solution (2.4). 



III. ROGUE WAVES IN THE 
D AVE Y- STEWART SON- 1 EQUATION 



As we have shown above, rational solutions (2.4) in 



Theorem 1 become rogue waves in the DSI equation when 
all parameters pi are required to be real. In this section, 
we analyze the dynamics of these rogue waves in detail. 



a = ai+ia2, b = b\ J rib 2l uj = uji + iuo 2l 9 = 9i J ri6 2l 
then 

2i(a 2 x + b 2 y + u 2 t + # 2 ) + 1" 



A(x,y,t) = V2 



W 



where 



W = (a 1 x + b 1 y + uj 1 t + 1 ) 2 + (a 2 x + b 2 y + uj 2 t + 2 ) 2 + /±. 



The solution Q can also be written down from (2.3) and 
the above /. 

This simplest rational solution has three distinctly dif- 
ferent dynamical behaviors depending on the parameter 
value of p\ . 

1. If p\ is not real, then it is easy to see that b/a is 
not real, hence b±/b 2 ^ a\/a 2 . In this case, along 
the [x(t),y{t)\ trajectory where 

a\x + b\y — —uj\t, a 2 x + b 2 y = —uj 2 t, 

(A, Q) are constants. In addition, at any given 
time, (A,Q) — > (a/2,1) when (x,y) goes to infin- 
ity. Thus the solution is a localized soliton moving 
on a constant background. 

2. If p\ < 0, i.e., pi is purely imaginary, then a, b and 
uj are all imaginary. In this case, the solution is a 
function of a 2 x + b 2 y + uo 2 t only, and is thus a line 
soliton moving on a constant background. 

3. If pi > 0, i.e., pi is real, then a, b are real but uj 
is imaginary. In this case, the solution is also a 
line wave, but it is not a moving line soliton any- 
more. As t ±oo, this line wave goes to a uniform 
constant background; in the intermediate times, it 
rises to a higher amplitude. Thus this is a line wave 
which "appears from nowhere and disappears with 
no trace" , hence it is a line rogue wave. 



A. Fundamental rogue waves 

The fundamental rogue waves in the DSI equation are 
obtained when one takes N = 1, n\ = 1 and p\ real 
in the rational solution (|2.4|), and en is a free complex 



parameter. After a shift of time and space coordinates, 
en can be eliminated and the fundamental rogue waves 
can be written as 



A(x,y,t) = V2 



1 



SiQt - 4 



1 + (k r x + k 2 y) 2 + 4QH 2 _ 



Q(x,y,t) = l-4ekl 



where 



fei = Pi - ePi 



2^2 



1 - (k ± x + k 2 yf 
[1 + (k ± x + k 2 y) 2 + 4m 2 } 2 



, (3.1) 



(3.2) 



: Pi + tPi 



ft = Pi + Pi 



This solution describes a line rogue wave with the line 
oriented in the {k\,k 2 ) direction of the (x,y) plane, 
thus the fundamental rogue waves in the DSI equation 
are line rogue waves. Along the line direction (with 
kiX-\-k 2 y fixed), the solution is a constant. As t — > ±00, 
the solution A uniformly approaches the constant back- 
ground a/2 everywhere in the (x, y) plane; but in the in- 
termediate times, \A\ reaches maximum amplitude 3 y/2 
(i.e., three times the background amplitude) on the line 
k\x + k 2 y = at time t = 0. This fundamental rogue 
wave is illustrated in Fig. [TJwith parameters e = 1 and 
Pl = 1.5. 

The above fundamental rogue waves in the DSI equa- 
tion are two-dimensional counterparts of the fundamen- 
tal (Peregrine) rogue waves in the NLS equation [4j [5] . 
Indeed, when we take e = 1 and p\ = 1 in the above 
fundamental rogue waves, we have k\ = 0, hence the so- 
lution A is independent of x and Q = 1. In this case, 
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FIG. 1: A fundamental rogue wave (3.1) in the DSI equation 
with parameters e = l,pi — 1.5. 



the DSI equation reduces to the NLS equation, and this 
fundamental rogue wave of the DSI equation reduces to 
the Peregrine rogue wave of the NLS equation. 

Non- fundamental rogue waves can be obtained from 
the i V-rational solutions of order (ni, ri2, • • • , tin) in Eq. 
(2.4) with real values of (pi,...,Pat) when N > 1, or 
ni > 1, or both. 



Below we consider two subclasses of 
these non- fundamental rogue waves. 



B. Multi- rogue waves 

One subclass of non-fundamental rogue waves is the 
multi-rogue waves, which are obtained when we take N > 
1, m = • • • = njv = 1 in the rational solution (2.4) with 



real values of (pi, . . . ,Pn)- These rogue waves describe 
the interaction of N individual fundamental rogue waves. 
When t — >• ±oo, the solution approaches the constant 
background uniformly in the entire (x, y) plane. In the 
intermediate times, TV line rogue waves arise from the 
constant background, interact with each other, and then 
disappear into the background again. In the far field 
of the (x, y) plane, the solution consists of N separate 
line rogue waves. However, in the near field where these 
line rogue waves intersect and interact, wavefronts of the 
solution are no longer lines, and interesting curvy wave 
patterns would appear. 

To demonstrate these multi-rogue- wave solutions, we 
first consider the N = 2 case. In this case, the / and g 
functions of the solutions can be obtained from (2.4) as 



m 

771 



m 21 rn 22 



(3.3) 



where 



m 



Pi+Pj 



m\ 3 



1 
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Pj 



(Pi +Pj) 



Pipj 



Pi+Pjj 1 (Pi+Pj) 2 

free real 



£j is given by Eq. (2.6), Pi,P2 are tree real param- 
eters, and en, C21 are tree complex parameters. The 
complex parameter en can be removed by a shift of the 
(x, y, t) axes, then this two-rogue-wave solution contains 
four non-trivial real parameters, namely, pi, P2, and the 
real and imaginary parts of C21. This solution for param- 
eters 



e = 1, pi = 1, p 2 = 1.5, en = c 2 i = 



(3.4) 



is shown in Fig. [2] It is seen that when these two line 
rogue waves arise from the constant background, the re- 
gion of their intersection acquires higher amplitude first 
(see t = — 1 panel). After these higher amplitudes in 
the intersection region fade, the line-rogue solutions in 
the far field then rise to higher amplitude (see t = 
panel). Interestingly, the wave pattern at t = features 
two curvy wavefronts which are well separated. These 
curvy wavefronts are caused by the interaction of the 
two fundamental (line) rogue waves. At large times, the 
solution goes back to the constant background again (see 
t = 5 panel). 

For larger N, these multi-rogue waves have quali- 
tatively similar behaviors, except that more line-rogue 
waves will arise and interact with each other, and more 
complicated wavefronts will form in the interaction re- 
gion. For example, with N = 3 and parameter choices 



e = 1, pi = 1, p 2 = 1.5, p 3 2, en 



C21 



C31 0, 
(3.5) 

the corresponding solution is shown in Fig. [3j As can be 
seen, the transient solution patterns become more intri- 
cate. 



C. Higher-order rogue waves 

Another subclass of non-fundamental rogue waves is 
the higher-order rogue waves, which are obtained when 
we take N = 1 and n\ > 1 in the rational solution (2.4) 



with a real value of p\. For ins tanc e, if ri\ = 2, we get 
second-order rogue waves from (|2.4|) as 
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FIG. 2: A two-r ogue - wave solution ( |3.3| in the DSI equation 
for parameters ( |3.4| ). Plotted is the p4f field. The constant- 
background value is y/2. 




FIG. 3: A t hree -rogue- wave solution in the DSI equation for 
parameters (3.5). Plotted is the \A\ field. 
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Pi +pi 



(3.7) 



where ^ is given by Eq. (2.6 ), pi is a free real parameter, 
and C12 is a free complex parameter. Here we have set 
en = in (2.4) by a shift of the (x,y,t) axes. Higher- 



order rogue waves with ri\ > 2 can be similarly obtained. 

An interesting phenomenon is that, unlike the multi- 
rogue waves discussed in the previous subsection, these 
higher-order rogue waves do not uniformly approach the 
constant background as t — >• ±00. Instead, only parts of 
their wave structures approach the constant background 
as t — » ±00, but the other parts move to the far distance 
as localized lumps with undiminished amplitude and in- 



creasing velocity as t —> ±00. To illustrate these behav- 
iors, we consider the above second-order rogue waves. 
For parameter values 



e = 1, pi = 1, C12 = 0, 



we get 



/ = - [x - {At 2 + y-y 2 



^2 1 

(y-- 2 f + - 4 



(3i 



(8t 2 + H ). 



g = f - 4it(x + y-y 2 - U 2 ) +x + y-y 2 - 12t 2 - -, 



and 



A(x,y,t) = V2^. 



(3.9) 



For this solution, A(x,y,t) = A(x,y,—t), thus 
\A(x,y,—t)\ = \A(x,y,i)\. This solution is displayed in 
Fig. [4} We see that when \t\ ^> 1, the solution is a 
localized lump sitting on the constant background \/2 
(see t = ±7 panels). The peak amplitude of the lump is 
attained at 



(x,y) = (4t 2 + 



thus this lump is accelerating rightward as \t\ increases 
(see t = — 7 and t = — 5 panels). The peak- amplitude 
value of the lump stays at 3^/2 and is unchanged though. 
As the lump accelerates rightward (with increasing |t|), 
its vertical size (along the y direction) remains the same, 
but its horizontal size (along the x direction) expands 
(see t = —7, —5, —1 panels) . When t — >• 0, this lump 
disappears. At the same time, a parabola-shaped rogue 
wave rises from the background (see t = — 1 and t = 
panels). At t = 0, this parabola is located at 

x = y-y 2 , 

where the rogue wave reaches peak amplitude 3y/2 (see 
t = panel). Visually one may describe the solution in 
Fig. [4] as an incoming lump being reflected back by the 
emergence of a parabola-shaped rogue wave. In addition, 
this lump decelerates as it comes afar and accelerates as 
it goes away. 

We have also examined the second-order rogue waves 
(3.6)- (3.7) for other parameter choices of (e,pi,ci2), and 
found that those solutions are qualitatively the same as 
the one in Fig. 4, except that those solution patterns 
may be stretched and skewed in the (x, y) plane. 



IV. SUMMARY 

In summary, we have derived general rogue waves in 
the DSI equation by the bilinear method, and our solu- 
tions are given in terms of determinants. We showed that 
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Lemma 1 Let m\j\ (p^ and be functions of xi, 
X2, X-i and x_2 satisfying the following differential and 
difference relations, 



(n) / (n) 



^ 2 m^ = ^ +1 Vf ) +^ n) ^ n " 1) , 



(n) 



(n-2) /(n+1) 



(n-1) ,(n+2) 



(n+1) (n) . (n) , (n+1) 



(n) 



(n+i/) 



(n) 



1,2,-1,-2). 



(A.l) 



FIG. 4: A second-order rogue wave solution (3.9) in the DSI 
equation with e = 1. 



Then the determinant, 



det 



K') . 



(A.2) 



the simplest (fundamental) rogue waves are line rogue 
waves which arise from the constant background with a 
line profile and then disappear into the constant back- 
ground again (see Fig. 1). We also showed that multi- 
rogue waves describe the interaction of several fundamen- 
tal rogue waves, and interesting curvy wave patterns ap- 
pear due to this interaction (see Figs. 2 and 3). However, 
higher-order rogue waves were found to show very novel 
features. Specifically, only parts of the wave structures 
in the higher-order rogue waves rise from the constant 
background and then retreat back to it, but the other 
parts of the waves come from the far distance as local- 
ized lumps, which interact with the transient rogue waves 
in the near field and then are reflected back and acceler- 
ate to the large distance again (see Fig. 4). These rogue 
wave solutions to the DSI equation generalize the rogue 
waves of the NLS equation into two spatial dimensions, 
and they could play a role in the physical understanding 
of rogue water waves in the ocean. 
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satisfies the bilinear equations, 

(D Xl D x _ 1 - 2)r n • r n = -2r n+ ir n _i, 
(D 2 X1 -D X2 )r n+1 -r n = 0, 
(D 2 X1 +D x _ 2 )r n+1 -r n = 0. 



(A.3) 



This lemma can be proved by the same method as for 
Lemma 3.1 in [12], thus its proof is omitted here. We 
note that by the variable transformation 

xx = -(x + y), x-! = ^(x-y), x 2 = - % -t, X- 2 = ^t 

(A.4) 

and the complex conjugate condition 

T~n T— n , (A- 5) 



the above bilinear equations (A.3) are reduced to the 



bilinear form (2.2) of the DSI equation for / = tq and 



g = t\. Therefore all we need is t o cho ose appropriate 



matrix elements rnffl which satisfy (lA.ll) and realize the 



conjugate condition (A.5) with (A.4). 



Proof of Theorem 1 It is easy to see that functions 
(f\ n \ ip^ an d mff defined by 



(n) 



(n) 
™>ij 



1 



(n) / (n) 7 

3 Pi + Qj Qj 



Appendix 



1 1 2 

£• = -^x-2 H X-i +PiXi J rp i x 2 , 

Vt Pi 



In this appendix we will prove Theorem 1 in section 
[n|by using the bilinear method. Firstly we present the 
following lemma. 



Vj 



x_2 H £-1 + qjXi - qjX 2 , 
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satisfy Eqs. (A.l). Here Pi : qj are arbitrary complex 
constants, and it is assumed that the lower boundary 
value of the integral in the above mf^ equation is zero. 
But these functions do not lead to rational solutions. 

To get rational solutions, we differentiate the above 
functions with respect to the parameters pi and qj. To 
obtain solution expressions in (a) of Theorem 1, we con- 
sider the following (p[ n \ an d m $ functions: 

V \ n) =Ai P ^, =B j (-q j )- n e^, (A.6) 



(-^)»e« 



Pi + qj qj 



(A.7) 



where Ai and Bj are differential operators of order rii 
and rtj with respect to pi and qj respectively, defined as 



Next we derive solution expressions in (b) of Theorem 



(n) (n) 

1. For this purpose, we consider the following (p\ , ip y J 



and m-^ functions: 

i^e* # ) = 4-(-«j)" n e % , (A.10) 



4? = p^^fdx,, 



(A.11) 



where and are differential operators of order 
and with respect to pi and ^ defined as 



4 = £ 



k=0 



(Hi - k)\ 



(M^, (A.12) 



k=0 



1=0 



(A.8) 

Qfc, (ij/ are arbitrary complex constants, and rii are arbi- 
trary positive integers. It is easy to see that these func- 
tions also satisfy Eqs. (A.l), thus r n = det (m^^ with 
(A.7) satisfies the bilinear equations (A.3). By using the 



operator relations 

( Pid Pi )p?e^ = p?e ?I ( Pi d Pi + & + n), 



4 = E 



1=0 



( n ._Z)|t<W ' 



(A.13) 



Cjfc, dji are arbitrary complex constants which are related 
to the constants dji in the solution expression (a) as 

Cik = c ik (ni - k)l, dji = dji(rij -1)1, 
and it is as sumed that the lower boundary value of the 



(A.10) as 



integral in (|A.ll|) is zero. Rewriting cp^ and in 



(?A)(-<&)-V» = (- qj )- n e^( qj d qj +r ] ' j -n), 



where 



2 


1 






Pi 


Pi 


2 


1 



3 <% 



X-2 x-i + qjX\ - 2q,X2, 

Qj 3 



the matrix element mffi in (|A.7|) becomes 



(n) 



(--) n e^^c, fc (^ i+ e-+n) r 
4/ 



A;=0 



i=0 



Taking parameter constraints 

qj=pj, dji=Cji, 



(A.9) 



(A.14) 



B J -(-«,-)- n e" J = Qf\-qj)- n e^, (A.15) 

where and are polynomials of degrees and 
rij in (x_2, xi, X2) respectively, and using the inte- 



gration by parts, mf^ in (|A.ll|) then turns into 



ni + n i ( n \i/ 



v=0 



Due to the gauge freedom of r n , we can see that 



(n) 



ni-\-rij 

E 



1 



Pi + Qj 



dl^Qf (A.16) 



gives the same solution. In order for this solution to 
satisfy the conjugate condition (A.5), we take 



qj=Pj, dji=Cji. 



(A.17) 



and using the variable transformation (A.4), we obtain 



jr (n) ( — n) 



thus the conjugate condition (A.5) is satisfied. Finally 
using the gauge freedom of r n , we obtain the rational so- 
lutions to the DSI equation as given in (a) of Theorem 1. 



Then in view of the variable transformation (A.4), we 
have 



]f=Pj~ n \ m£ ) =m<r n) , f n = r- n , (A.18) 



hence r n with the above matrix elements (A.16) satisfies 



the bilinear DSI equation (2.2). 



8 



Lastly we derive the explicit expression of . From 



the definition of i^ (n) in dA.14l, we have 



£ 

k=0 



(n t -k)V 



For £ = p v x vi using the functional identity 

^=-2,-1,1,2 

(see Appendix of [12]) 



we get 
1 



e K ^F(p) = F(e K p) 1 



where £ (n) (p) = (&(p) + rc,6(p), • ■ • ,&(p) + 4m, • • • ) 
and £fc(p) = y k p^XjJ_k \. By comparing the coefficient 
of order K k in Eq. (A.20), we obtain 



1 (Pdp)\n 



f-p n e* = S fc (£ (n) (p)). 



p n e^ fc! 



Substituting p = ^ a nd £ = ^ into this equation and 
inserting it into (A. 19), and recalling the variable trans- 



J> (n) as 



formation (A.4), we then find the explicit expression for 



rii 



00 k ( \ 

exp (nn + — v k p v x \ where £ w (pi) is as given in Th eorem 1. Inserting this 

V ^ k\ ^ V pM anc i QW f rom (|A.18| into Ifil t.ViP solution py- 



fc=l ^ 

CO 

V> fc S fe (£ (n) (p)), (A.20) 

k=0 



pression in (b) of Theorem 1 is t 
the proof of Theorem 1. 



A.16), the solution ex- 



ren derived. This ends 
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